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We report on an outburst of the high mass X-ray binary 4U 01 15+63 with a pulse period of 
3.6s in spring 2008 as observed with INTEGRAL and RXTE. By analyzing the lightcurves we 
derive an updated orbital- and pulse period ephemeris of the neutron star. We also study the 
pulse profile variations as a function of time and energy as well as the variability of the spectral 
parameters. We find clear evidence for at least three cyclotron line features. In agreement with 
previous observations of 4U 01 15+63, we detect an anti-correlation between the luminosity and 
the fundamental cyclotron line energy. 
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Figure 1: Left: RXTE ASM lightcurve of the outburst in spring 2008 with a binning of 0.5 days together 
with the observation times of RXTE (green) and INTEGRAL (blue). Right: Doppler shifted pulse periods of 
4U 01 15+63 versus the orbital phase together with the best fit (<p = 0 where pulse period is at maximum). 

1. Introduction 

4U 01 15+63 is a transient X-ray pulsar with a Be star optical companion. Unlike ordinary B 
stars, Be stars host a circumstellar disc, resulting from fast rotation, magnetic loops, or non-radial 
pulsations. Due to the highly eccentric orbit of the system, material originating from a Be star’s 
disc can be accreted onto the neutron star during the periastron passage. One of the main processes 
leading to the observed hard X-ray emission is the Comptonisation of seed photons, e.g., from 
blackbody radiation of a hot spot on the neutron star surface, in the hot plasma of the accretion 
column above the magnetic poles. 

On the surface of neutron stars magnetic fields of the order of 10 12 G can occur. As a conse- 
quence, the X-ray spectra of pulsars may contain cyclotron resonance scattering features (CRSF) 
at a fundamental line energy of 
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as well as at the associated harmonic energies, where z g is the gravitational redshift. So far CRSFs 
have been detected in about 20 X-ray pulsars. 

4U 01 15+63 is one of the pulsars for which CRSFs have been studied in great detail (see, e.g., 
[4, 5,8,9, 11, 13, 14]). In previous outbursts, CRSFs have been detected up to the fourth harmonic 
[11]. This high number of detected CRSFs in 4U 0115+63 makes this system an outstanding 
laboratory to study the physics of cyclotron lines in X-ray pulsars. 

The overall lightcurve of the outburst we report on is shown in Fig. 1 (left). The event took 
place in spring 2008, with the peak intensity amounting up to ~280mCrab. The outburst was 
observ ed in 21 individual RXTE- and 8 INTEGRAL pointings. 
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Figure 2: Left: Pulse period ephemeris based on the binary-corrected pulsar lightcurves derived by phase 
connection (top panel) and the phase shift between the measured pulse profiles and the pulse period 
ephemeris (bottom panel). Right: P versus the luminosity. The colours red, purple, and blue correspond 
to the brightening phase, the outburst’s maximum, and the dimming phase, respectively. Both figures are 
preliminary. 

2. Orbital Ephemeris 

Due to the orbital motion of the neutron star, its pulse period is Doppler shifted periodically 
(see Fig. 1, right). These shifted periods, calculated by epoch folding the PCA-lightcurves, can be 
used to refine the orbital ephemeris of the binary system. We used the orbital parameters from the 
literature [1] (a x sini = 140.13(8) lt-sec, e = 0.3402(2) and a> = 47.66(9)°) except the epoch of 
periastron passage. To, which we determined from the fit. Assuming that the orbital period P or b 
remains constant, the difference between our To and that from [1] has to be an integer times P 0 rb- 
This leads to P or b = 24.3 16171^^^1 d. Recently, evidence for a non zero time derivation of the 
longitude of periastron do has been found [10]. The study of how this evolution affects our results 
is still work in progress. 


3. Pulse Period Ephemeris 

Using the new orbital ephemeris, we corrected the PCA-lightcurves for binary motion. The 
remaining variations in the pulse period P p are due to transfer of angular momentum of the accreted 
matter onto the neutron star. By epoch folding these corrected lightcurves and performing phase 
connection (see, e.g., [12]), we calculated the pulse period ephemeris displayed in Fig. 2 (left). The 
respective parameters are summarized in Table 1, while the pulse period at a certain time t can be 
calculated as 

P(t) = Pp(7b,p) + Pp(t - 7o, p) + ~P p (t - 7b.p) 2 . (3.1) 

As the first and the second time derivative of the pulse period, P p and P p , are expected to change 
significantly during the outburst, we analyzed four different time intervals seperately. As displayed 
in Fig. 2 (left), the evolution of the pulse period is not continuous over the outburst. Especially 
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Figure 3: Pulse profiles of 4U 0115+63 for several energy bands and times. 


the strong spin down-"jump” between the first and the second time interval is still not understood. 
Unfortunately there are only three data points available during this epoch. Because the times when 
the respective pulse profiles have been measured are too far away from the neighboring intervals, 
P p and P p at these three points are expected to be different than during the intervals with a well 
constrained ephemeris. Therefore we had to omit these three data points from our analysis. Despite 
this enigmatic behaviour, the lower panel of Fig. 2 (left) shows that the phase shifts between the 
observed pulse profiles and the pulse period ephemeris (see Table 1) are small enough such that 
this pulse period ephemeris is suitable for phase resolved spectroscopy. 

Fig. 2 (right) shows that the 3-50 keV luminosity and P p are correlated for the majority of data 
points. This is in agreement with the classical interpretation (see, e.g., [3], and references therein), 
where for higher luminosities a higher mass transfer rate and therefore a higher transfer of angular 
momentum is expected. 
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Table 1: Preliminary parameters of the pulse period ephemeris. The epochs are marked as follows: (i) 
MJD < 54552.0, (ii) 54557.0 < MJD < 54564.5, (iii) 54564.5 < MJD < 54575.0, and (iv) MJD > 54575.0 


epoch 

phase 

Fo.p (MJD) 

Fp(To.p)[Sj 

/> p [ss l ] 

Pp[ss- 2 ] 

H4 

(i) 

brightening 

54550.055284 

3.614200 

-1.79 x l(r 10 

-4.96 x 10~ 16 


(ii) 

maximum 

54558.948604 

3.614260 

-6.24 x 10- 10 

+ 1.75 x 10~ 15 

Us 

(iii) 

dimming 

54570.138089 

3.614120 

-3.43 x lO^ 11 

+6.22 x 10- 17 

(iv) 

dimming 

54578.350020 

3.614125 

±0 (fix) 

±0 (fix) 



4. Pulse Profile Variations 

Pulse profiles for different energy bands and epochs are shown in Fig. 3. In agreement with 
previous work [ 1 ] the secondary peak is more prominent at lower energies, while at higher energies 
it disappears almost completely. These two peaks in the pulse profile are thought to originate from 
the two hot spots at the magnetic poles of the neutron star. An offset dipole field could explain 
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Figure 4: 3-50 keV luminosity versus fundamental cyclotron energies. The small box provides a close-up 
view for high luminosities. For L\ > 9 x 10 37 ergs _1 the fundamental line energy stays nearly constant. 

the different behaviour of the two pulse peaks, which has been observed in other transient X-ray 
sources as well (e.g., GRO J1008-57: [6], or 4U 1909+07: [2]). 


5, Spectral fitting 


We modelled PCA- and HEXTE data taken during the maximum of the outburst using the 
NPEX model (see., e.g. [7], and references therein), which can be written as 


NPEX(£) = {A\E~ a ' +A 2 E +a -) x exp 



(5.1) 


where Ain are the normalization constants, a\j 2 are the photon indices of the positive and negative 
powerlaws, and kT represents the cutoff energy. We included three cyclotron absorption features 
at ~ 1 1 , 22, and 33 keV using the model CY CLABS, which is given as 


CYCLABS (E) = exp 


( r„(W n E/E„) 2 \ 

V (E-Erf + Wi)’ 


(5.2) 


where E n is the nth resonance energy, W n the corresponding width, and x„ the optical depth. Fur- 
thermore we included a Gaussian iron emission feature at 6.4 keV and achieved an acceptable fit 
with Jred ~ 1.1. Using these CRSF- and continuum parameters as a starting point and keeping 
(7.2 —2 , Eo = 2E\ , and Wo = Wj fixed, vve modelled all available spectra of the outburst. For all 
observations we obtained between 1 and 2, but mostly less than 1.5. The photon index <*[ 
turned out to be constant within the uncertainties during the outburst and was therefore fixed to 
1.2 to constrain the other parameters better. We determined values of 6.0±0.2keY r for the cutoff 
energy kT in the beginning, and 4.6±0.2keV at the end of the outburst. While we found three 
cyclotron line features in all spectra until MJD=54568, after this time only the fundamental line 
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was detected. An anti -correlation between the fundamental cyclotron line energy and the 3-50 keV 
luminosity has been found for a previous outburst in 1999 [9]. As displayed in Fig. 4, this result is 
consistent with our new analysis. The hysteresis-like shape of the curve might be due to a change 
in the accretion geometry during the maximum of the outburst. This might be also the reason for 
the jump in the pulse period ephemeris of the neutron star detected during the maximum of the 
outburst. 

6. Outlook 

As forthcoming work we will check the consistency of the predicted changes in the periastron 
longitude [10] with our results. Furthermore we will analyze the phase resolved behavior of the 
spectral parameters, based on our updated pulse period ephemeris. For this purpose we will also 
include the INTEGRAL data into our spectral analysis. 
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